MZ-TH/10-03 



Ghost wave-function renormalization in 
Asymptotically Safe Quantum Gravity 

Kai Groh and Frank Saueressig 

Institute of Physics, University of Mainz 
Staudingerweg 7, D-55099 Mainz, Germany 

kgrohOthep . physik .uni-mainz . de 
saueressigSthep . physik . uni-mainz . de 

Abstract 

Motivated by Weinberg's asymptotic safety scenario, we investigate the gravitational renormal- 
ization group flow in the Einstein-Hilbert truncation supplemented by the wave-function renor- 
malization of the ghost fields. The latter induces non-trivial corrections to the /3-functions for 
Newton's constant and the cosmological constant. The resulting ghost-improved phase diagram 
is investigated in detail. In particular, we find a non-trivial ultraviolet fixed point, in agreement 
with the asymptotic safety conjecture which also survives in the presence of extra dimensions. 
In four dimensions the ghost anomalous dimension at the fixed point is t]* = —1.8, supporting 
space-time being effectively two-dimensional at short distances. 



1 Introduction 



Constructing a consistent and predictive quantum theory for gravity is one of the prime chal- 
lenges of theoretical high energy physics today. One proposal in this direction, which recently 
received a lot of attention, is Weinberg's asymptotic safety scenario [1, 2, 3], see [4, 5, 6, 7, 8, 9] for 
reviews. This scenario adopts Wilson's modern viewpoint on renormalization [10], assessing that 
gravity has a fundamental description within the framework of non-perturbatively renormalizable 
quantum field theories. The key ingredient underlying this idea is a conjectured non-trivial (or 
non-Gaussian) fixed point of the gravitational renormalization group (RG) flow. For RG trajec- 
tories attracted to it at high energies, the fixed point ensures that the dimensionless coupling 
constants remain finite, so that physical quantities are safe from unphysical UV divergences. 
These trajectories, also including the one describing our world [11], span the UV critical surface 
Sjjy of the fixed point. The precise position of "our RG trajectory" within this surface is deter- 
mined by free parameters, which have to be fixed by experiment. For a finite-dimensional 5uv 
this construction is then as predictive as a standard, perturbatively renormalizable, quantum field 
theory. Elucidating the fixed point structure of the gravitational RG fiow and understanding the 
properties of the corresponding UV critical surfaces is a central aspect of the asymptotic safety 
program to date. 

An important technical tool in this program is the functional renormalization group equation 
(FRGE) for gravity [12]. Formulated in terms of the Wetterich equation [13], the FRGE describes 
the dependence of the effective average action on the coarse graining scale k 



Here, F^ ' denotes the second functional derivative of the effective average action with respect to 
the fields of the theory and STr is a generalized functional trace which includes a minus sign for 
ghosts and fermions. Furthermore, TZk is a matrix valued IR cutoff, which provides a fc-dependent 
mass-term for fluctuations with momenta < k'^. The interplay between the full regularized 
propagator ^F^,^^ -|- TZ^j and dfTZk ensures that the STr receives contributions from a small 
p^-interval around k"^ only, rendering the trace-contribution finite.^ 

While (1.1) constitutes a formally exact equation, it comes with the drawback that it cannot 
be solved exactly. A widely used approximation scheme for obtaining non-perturbative infor- 
mation from the FRGE consists of making an ansatz (truncation) for T/. which retains a finite 
number of A;-dependent coupling constants only. Projecting the the resulting RG flow onto the 
truncation subspace allows to read off the /3-functions for the running couplings as the coef- 
ficients of the interaction monomials retained by the ansatz without resorting to perturbation 
theory. Employing truncations, an important point is to establish the robustness of the emergent 
physical picture. Within a given ansatz, its reliability may be tested by studying the dependence 
of physical quantities on the shape of the unphysical IR regulator TZ^. A second, more laborious, 

^The FRGE captures the RG-dependence of the effective action. The relation between the effective action at a 
fixed point F* and the corresponding fundamental action S has recently been discussed in [14]. 
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route consists in extending the truncation ansatz, thereby showing that all physical results, as, 
e.g., the fixed points of the RG flow, are robust under the extension. 

For gravity, these ideas have been implemented systematically, so far focusing on the gravita- 
tional sector of Fjfc mostly. In this class the truncation ansatz is spanned by diffeomorphism invari- 
ant operators build from the physical metric g, e.g., = (167rGjfc)~^ / d'^Xy/g{—R + 2Afc) + 

The most studied case, the so-called Einstein-Hilbert truncation, encompasses a scale-dependent 
Newton's constant and cosmological constant Aj^., and has been analyzed in a number of 
works [15, 16, 17, 18, 19, 20, 21]. Subsequently, this ansatz has been refined by including higher- 
derivative i?^-intcractions [22, 23, 24], higher order polynomials in R up to within the frame- 
work of /(i?)-gTavity [25, 26, 27], non-local operators [28, 26], and lately also the Weyl-squared 
interactions [29, 30, 31] capturing the characteristic features of higher-derivative gravity. All 
these computations have identified a NGFP of the gravitational RG flow, providing substantial 
evidence for the asymptotic safety scenario. Moreover, the f{R)- and C^-results point at the di- 
mension of the associated UV critical surface being finite, possibly even as low as three. Notably, 
the essential features of this picture already emerge from the structurally significantly simpler 
flow equations obtained within the conformally reduced gravity framework [32, 33, 34]. 

In this work we go beyond the gravitational approximation of F^, including quantum effects 
from the ghost sector. More specifically, we will augment the Einstein-Hilbert truncation by a non- 
trivial wave-function renormalization of the ghost fields (for a similar computation in Yang-Mills 
theory see [35]). Our main motivation for this ghost-improvement originates from the analogy to 
QCD where this coupling plays an essential role for the IR physics of theory [36, 37, 38, 39]. While 
it is clear, that there is also a non-trivial interplay between the ghosts and the gravitational /3- 
functions in gravity, explicit computations are not available yet. Here we close this gap, computing 
the anomalous dimension of the ghost fields r)c and its effect on the running of Newton's constant 
and the cosmological constant. The explicit computation is based on a new perturbative heat- 
kernel technique developed by Anselmi and Benini [40] which allows the systematic expansion 
of the flow equation in the presence of background ghost fields in a curved background.^ This 
continues the exploration of quantum gravity effects in the ghost sector [42]. 

The rest of the paper is organized as follows. In Section 2 we derive the /9-functions governing 
the RG dependence of Newton's constant, the cosmological constant, and the wave-function 
renormalization in the ghost sector, using a perturbative heat-kernel technique for non-minimal 
differential operators. The ghost-improved fixed point structure and phase diagram is analyzed in 
Section 3 and we comment on our findings in Section 4. Some technical details on the heat-kernel 
techniques and threshold functions employed in this paper are relegated to Appendix A, while the 
functions determining the gravitational and ghost anomalous dimensions are defined in Appendix 
B. 

^In a companion paper [41], a similar analysis will be carried out using a spectrally adjusted cutoff and a 
flat-space projection technique. We thank A. Eichhorn and H. Gies for informing us on their upcoming work. 
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2 The ghost-improved Einstein-Hilbert truncation 



The main purpose of this paper is the investigation of the gravitational RG flow, including the 
quantum effects captured by the wave- function renormalization in the ghost-sector. Our trun- 
cation ansatz, which will be called the "ghost-improved Einstein-Hilbert truncation" , encompasses 
three scale-dependent coupling constants: Newton's constant Gk, the cosmological constant A^, 
and the power-counting marginal multiplying the ghost-kinetic term. In this section, we will 
derive the non-perturbative /3-functions capturing the RG dependence of these couplings. 

2.1 The truncation ansatz 

Our ansatz for the effective average action is of the general form 

r^b, c, c- g, c, c] = rr^b] + rf b; g] + rf[g, c, c- g, c, c] . (2.1) 

Besides on the physical metric g and the classical ghost fields C,C, it also depends on the 
corresponding background fields g and c, c. They are related by 

Qui' — 9iJ.i> + h^ii/ , C/j, = Cn + f/j, , Cn = C/j, + fn , (2-2) 

where h/^^, and denote the expectation value of the quantum fluctuations around the back- 

ground, which are not necessarily small. In the gravitational approximation the computations 
are simplified by setting the background ghost fields to zero. This, however, does not allow to 
keep track of the ghost-kinetic term, so that we work with a non-trivial ghost background in the 
following. 

The gravitational part r|'^'^^[g'] is build from the physical metric g and taken to be of the 
Einstein-Hilbert form 

rr^'b] = 2 Jd'^x^ {-R + 2Afc) , (2.3) 

where = (327rGo)~^ with Go a fixed reference scale and Z^ denotes the wave-function renor- 
malization for the graviton. r|'^'^^[5r] is supplemented by the gauge-fixing term ^^[g', §]■ Employing 
the harmonic gauge, the latter reads 

rf [h; g] = \Z^ j d^xV^ g^'F^^Fu , = V2«:(^^ ^ " ^^/.^^''^a/?) • (2-4) 

The resulting Faddeev-Popov determinant is captured by the ghost term 

rf [5, C, C- g, c, c] = -^Zl J d^x^ A^^ , (2.5) 

with 

A^M^ = r^r^DxigpuDa + ga.Dp) - gP'^gi'^Dxiga.Dp) , (2.6) 

containing the wave-function renormalization of the ghosts Z^. The gauge-choice (2.4) has the 
main virtue, that it allows for a straightforward comparison to earlier results [15] obtained in the 
Einstein-Hilbert truncation without ghost-improvement, fixing Z^ = 1. 
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Before entering into the explicit computation of the ^-functions arising from the ansatz (2.1), it 
is useful to first consider the left-hand-side of the flow equation and identify the interaction mono- 
mials whose coefficients encode the running of our coupling constants. Taking the Sj-derivative 
of our ansatz (2.1) and setting the fluctuation fields to zero afterwards yields 

dtVj, = 2k^ I d^xVS HdtZ^)R + ^dtiZ^Ak)] - V2{dtZ',) J d'^Xy/^ c^D^c^ + • • • . (2.7) 

Thus it suffices to extract the Einstein-Hilbert monomials and the ghost kinetic term from the 
right-hand-side of the flow equation. Comparing the coefficients multiplying these interactions 
then allows us to read off the desired /3-functions for Z^, and Aj.. 

2.2 Quadratic forms and the inverse Hessian {t'-^^ +nk)-^ 

Upon specifying our truncation ansatz, we proceed by computing the second variation of Fj. with 
respect to the fluctuations (2.2). To facilitate the subsequent steps, it is convenient to decompose 
the metric fluctuations into their traceless and trace part 

The forms quadratic in the fluctuations can be simplified further by first identifying g = g,C = 
c,C = c and specifying a particular class of backgrounds. In our case, this class has to be 
general enough to distinguish the interaction monomials in the Einstein-Hilbert term and the 
ghost-kinetic term. This can be accomplished by choosing the background metric g as the one of 
a d-dimensional sphere, implying 

Moreover, the background ghost field can be taken transversal 

D^c>' = 0, D^c^ = 0, (2.10) 

which suffices to keep track of the ghost kinetic term. In the sequel, we will resort to this choice 
of background to simplify all expressions. For later reference, we also introduce the unit on the 
space of traceless symmetric tensors (2T) and vectors (1) 

l2T = i (5^5^ + 6;6P) - iy^,/- , li = ,5;; . (2.11) 

Given these preliminaries, we now expand the ansatz (2.1) around the background (2.9) and 
(2.10), retaining the pieces quadratic in the fluctuations only. For F|'^'^^[5] -|- r^[h,g] the result 
has already been given in [12] 

rr^[^; 9] = ^^Z^ Jd^xV^ { V [A - 2Afe + CtR] h'^" -^h[A- 2Afe + CgR] h} (2.12) 
with 
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and A = — 5^ denoting the covariant Laplacian constructed from the background metric. Owed 
to the non-trivial background ghost field, the analogous computation for r|'' is slightly more 
involved. Thus we first give the intermediate results obtained from expanding (2.6) contracted 
with a dummy vector in h^n 

M^, C + ^ R] C , 



(2.16) 



Here the bar indicates that we expand around g = g- Remarkably, the last line is independent of 
the choice of background or gauge. Based on these results, we obtain the quadratic form in the 
ghost sector 

(2.15) 

Here the two lines are equivalent up to surface terms and define the Grassmann-valued operators 
A, A, Q, Q, and their adjoints, respectively. Their explicit expressions read 

i« = 1 [D^c'^ + i2°^c^ + D^c'^D" + (2 - ^)D'^c''D^ + (1 - ^yD'^D^] , 

^« =i [5,c"L»" + 5"c"5, - (1 - f )c'^5<,5"] , 

and 

- ^gf.uiD^c'' + R'^ca + D'^c'^D, + iD'^c^D" + c'^D'^D^), 

+ ^r^iD'^CaD, + DaCaD-), (2-17) 

- \g^, [-CaD'^D'' + D'^c'^D^ + D'^d'D^] , 
Q^i^'^ = - \D„(f^D'' + D'^c''^ + D'^h^D'' + R^^^^c^ + D^^'c'^Wa 

+ [D^CaD^ + D^Ca + DaC^D'' + Rc.aC"] , 

where the covariant derivatives to the left of the ghost fields act on c or c only and {^u) = 
^{fiv + lyfi) denotes symmetrization with unit strength. The last lines in the Q-expressions 
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ensure that the operators are traceless in the 2T-indices fj,^. The quadratic forms (2.12) and 
(2.15) provide the key ingredients for constructing the IR cutoff TZf- and the inverse Hessian 
(rf ^ + 7^fe)-Mn the sequel. 

The construction of the IR-cutoff follows the prescription in [12]. I.e., at the level of the path 
integral, the gauge-fixed action is supplemented by an IR regulator of the form 

AkS = J (fxV^ { i V J/ Kp + \h TZkM h + T T^kjf U) ■ (2-18) 

The matrix 7?.^ is designed in such a way, that it adds a fc-dependent mass-term to the Laplacians 
appearing in the kinetic terms, 

A ^ Pfe = A + , (2.19) 

where Rk = k'^R^^\A/k'^) is the scalar part of the IR regulator and R^^ a dimensionless shape- 
function interpolating monotonically between R^\(}) = 1 and lim^-^oo ^^'^^ (-2) = 0. Comparing 
(2.18) to the quadratic forms (2.12) and (2.15), this prescription fixes 



7^ 



k,hh 



n 



a/2 li ZtRk 



(2.20) 



Observe that TZ^ inherits a non-trivial /c-dcpcndence via the wave-function renormalization fac- 
tors. Notably, it is this feature which is partially responsible for the non-perturbative character 
of the computation. 

(2) (2) 
The next step consists in constructing the Hessian and the inverse of F^, + TZk ■ Here we 

follow the conventions [37, 35], using a skew-symmetric metric to couple anti-commuting fields 

(see Appendix A of [37] for details).^ In terms of the multiplets 



$ = {V. h, r , fa}, ^ = { v> ^> -fa, n , 

the Hessian is given by the following matrix 



^(2) 



where all variations act from the left. Substituting the quadratic forms (2.12) and (2.15), 



(2.21) 



(2.22) 



(rf + 



K Q 
Q M 



assumes block form with entries 



K=K^ diag [{Pk - 2Afc + CtR)12T , -^{Pk - 2Afe + CgR)] 
M=V2ZI {Pk - iR) diag [ li , li ] , 



and 



Q = V2Z, 



Q = ^Z. 



(2.23) 



(2.24) 



(2.25) 



'One can explicitly verify that the resulting FRGE is equivalent to [18]. 
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The inverse of (2.23) is then found via the general inversion formula for 2 x 2-block matrices 



" K 


Q ' 


-1 


Q 


M 





-M-^Q (^K - QM-^qy 



-K-^Q [M -QK-^Q 

[m-qk-^qY^ 



-1 



(2.26) 



Expanding the inverse up to second order in the background-ghost fields, and taking into account 
the minus sign originating from the ghost sector of the supertrace, the right-hand-side of the flow 
equation becomes 



dt^k — ^Trgr 



[K- 



with 7^f 



- iTVgh [(M-^ + M-^QK-^QM-^)dtnf 
'■= + Sq + Si + Q2T + Go + Qi-\ , 



(2.27) 



diag[7^^ ^^,7?.fc^/j/j], TZf. = diaglTZj^ jf^TZj^ ^j], and terms not contributing to the 
truncation indicated by the dots. As already anticipated in the last line, the full trace decomposes 
into operator traces on the space of traceless symmetric tensors (2T), scalars (0) and vectors (1). 
Substituting the explicit expressions for the block matrices, the traces Si, which by definition, do 
not include Q and Q and are thus independent of the background ghost fields, become 

I2T 



52T = 5Tr2T 

So = iTro 
S, = - Tn 



[Z^iPk-2Ak + CTR) 



1 



Zj^iPk - 2Ak + CsR) 
li 



(2.28) 



zm - iR) 

They give rise to the /3-functions for Newton's constant and the cosmological constant. The f5- 
function for is captured by the terms of second order in the background- ghost fields. Neglecting 
the curvature terms and making use of the cyclicity of the trace, these are found as 



G2T 
Go 



V2k2(Z,^)2 



Tr2T 



2d 



d-2 V2av2(Z^^^2 



TVn 



dtjZj'Rk) 
[Pk - 2A]2 

dt{Z^Rk 



a ji Q per 

-He 



a Q per 
Pk 



Gi=- 



'k 

rTri 



[Pk - 2A]2 



Pk 



p2 
^k 



Qa 



1 



Pfc-2A 



V2«2Zf 



0, 



Pk 



(2.29) 



Pfc-2A 



Q 



1x1/ 



Notably, the insertions including the background ghosts, like, e.g., the QP~^Q and —QP~^Q 
in G2T, always appear in pairs. Each part thereby gives exactly the same contribution to the 
Gi, consistent with the symmetry factors of the underlying Feynman-diagrams. Evaluating both 
parts individually provides a highly non-trivial crosscheck when computing the coefficients of the 
ghost-kinetic term in the next subsection. 
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2.3 The ghost-improved /3-functions 

Wc arc now in the position to construct the /3-functions for , Z^, and Ak- In this context, it is 
useful to introduce the anomalous dimensions of Newton's constant and the ghost wave-function 
renormalization 

VN = -dt ln(Zf ) , r/e = -Ot In(Z^) , (2.30) 



together with the dimensionless couplings 

gk = Gkk''-^ = iZ^r'Gok''-\ 



Afe = Akk' 



(2.31) 



As explained in Section 2.1, the /5-functions for Newton's and the cosmological constant are 
encoded at zeroth order in the background-ghost fields. The desired interaction monomials are 
generated by the traces Si, eq. (2.28), which are evaluated straightforwardly by applying the 
early-time heat-kernel expansion detailed in Appendix A.l followed by inserting the identity 
(A. 15). Equating the resulting coefficients with (2.7) provides the flow equations for A^ and 
Gk = (327rZf 



1 



167rGfc 



k'i f d{d+l) 
(47r)<i/2 \ 2 

\ _ k''-^ ( d{d+l) 
J (47r)rf/2 \ 12 



^d/2 

a,o 



??^l>iy°(-2Afc) 

1^ ^1.0 



2d 



^i/2(0) 



^d/2-i(-2Afe)-iT?iv$,;2-i(-2Afe) 



d{d-l) 
2 



2,0 
d/2^ 



-2Xk) 



1 ^2,0 / 
2VN^d/2i 



-2Xk) 



1,0 
d/2- 



.i(0) 



2'lc^d/2- 



i(0) 



,^d/2(0) 



l^c*5/° (0)] } 



(2.32) 



In the limit Z^ = l,ric = 0, this result agrees with earlier computations [12, 15]. The terms 
proportional to r)c are novel and capture the backreaction of the quantum effects in the ghost 
sector on the running of the gravitational coupling constants. 

The final step is the computation of r]c- This requires extracting the background-ghost kinetic 
term from (2.29). Unfortunately, the differential operators entering into these traces arc not mini- 
mal, so that the early-time hcat-kcrncl expansion is no longer applicable and a more sophisticated 
method is needed. Here we follow Anselmi and Benini [40], which allows to compute operator 
traces of Laplace-type operators which include arbitrary (non-Laplace type) "perturbations" by 
considering the heat-kernel expansion at non-coincident points. Adapting these techniques to 
the case where the perturbation is build from the background-ghost fields, this method becomes 
applicable to Qi. (See Appendix A.2 for more technical details.) Retaining the ghost-kinetic term 
only, a tedious but straightforward computation yields 



G2T = 
00 = 
Gi = 



V2 Zl 
(47r)'i/2 K2(Zf )2 

V2 Zl 

1 1 

(47r)d/2 V2k2Z^ 



4d'^-d-8, 



4, Qd/2+l[Wl']-^Q,/2+2m' 

^Qd/2+i[Wf]-'^Qd/2+2[Wi']\ I d''xVSc^D%, 

2d2_5d_2^_, rwci I ^r^_ rwcil / ^d^ /^^Mn2, 



(2.33) 



'^2^Qd/2+i[Wt] +dQa/2+2m] / d'^xV^g-c^^D'c^ . 
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Here the Q-functionals are defined in (A.4) and 



TT.JV ^ dtjZ^Rk) j_ 
""l - (Pfc-2A)^ Tk 



""2 - (Pfc-2A)^ dx^;^ 



^5 



— ^2 



dxVk{x)-2A 



(2.34) 



a;=A 



Substituting (A. 16) and equating the result with the ghost kinetic term in (2.7), we finaUy find 



dfZ, 



(47r)'^/2- 



^2Ai 



+ Cgh (^i/2+i(-2Afe) - ^%|.;;,^,(-2Afe)) 

+ C/(W - 77c) (^'rf/2+2(-2Afc) + irf/2+2(-2Afe) 



(2.35) 



where 



4^2 - 9d - 2 



2^2 - 5d - 2 



Here the ^^/2+2 ^'^'^ ^d/2+2"*®^™^ originating from (5^/2+2 [^^2] drop out, due to the cancellation 
of the corresponding coefficients. 

The /3-functions are then obtained by solving eqs. (2.32) and (2.35) for dtXk, dtGh ^^nd 7]^- 
This yields 

dt^k = Px , dtgk = Pg = {d-2 + r]N)gk, 

with 

(2 - r/jv)Afc + hki^T^f-'^''' 



(2.37) 



/3a 



2d{d + 1)$^;,° (-2Afe) - ^d^dU^) - d{d + l)m^'^/2i-'^h) + 4d%$d/2(0) 
and the anomalous dimensions 



(2.38) 



Vn 





C3(A)C4(A)-5i(A)C2(A)" 




1-5 


52(A) + C2 (A) 


+ 52 


i?2(A)C2(A)-Ci(A)C3(A)" 



5C4(A) + 5' 


Si(A)Ci(A)-S2(A)C4(A)" 




1-5 


S2(A) + C2(A)" 


+ 52 


's2(A)C2(A)-Ci(A)C3(A)" 



(2.39) 



The explicit form of the functions i?j(A) and Cj(A) is given in Appendix B. The eqs. (2.38) and 
(2.39) are the desired /3-functions governing the RG dependence of gk,Xk, and rjc and constitute 
the central result of this section. 

Some comments are now in order. The inclusion of the wave-function renormalization for 
the ghosts gives non-trivial contributions to the /3-functions for g, A. These encompass the terms 
proportional to rjc in Px and the qualitatively new g'2-terms in t/at- The results obtained within the 
standard Einstein-Hilbert truncation [12, 7] are recovered by setting Ci = 0. Comparing powers 
of g, the leading contributions from the ghost sector are suppressed by one power of g, relative to 
the leading Einstein-Hilbert terms. Thus in the classical regime, g <^ 1, the ghost-improvement 
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may be neglected. In the quantum regime close to the NGFP where 5 ~ 1, however, we expect 
that these corrections become important. Investigating the influence of these new terms on the 
gravitational RG flow is then the subject of the next section. 

3 Properties of the flow equation 

After deriving the /3-functions (2.37) and (2.39), we now proceed by studying their properties, 
mostly resorting to numerical methods. In this context, it is useful to observe that Z£ enters into 
(3\ via ?7c only and is in turn completely determined by gk,)^k- Thus, substituting the explicit 
formula for tJc into Px, the running of decouples and allows to analyze the gravitational RG 
flow in the two-dimensional g, A-subsystem. Once a RG trajectory for g^, is found, it can be 
plugged back into r)c to obtain the running of the ghost anomalous dimension. We will now exploit 
this decoupling and first discuss the fixed point structure of the ghost-improved Einstein-Hilbert 
truncation for c? = 4 in Section 3.1, before focusing on the phase portrait and the fixed point 
structure including extra-dimensions in Sections 3.2 and 3.3, respectively. 

3.1 Fixed points of the four-dimensional /3-functions 

As highlighted in the introduction, the crucial ingredient of the asymptotic safety scenario is 
the fixed point structure of the gravitational /3-functions. Thus, we start our investigation by 
looking for fixed points g*,\* where /3g = (5x = simultaneously. In the vicinity of such a 
point, the linearized /3-functions are given by dtgi = Bij{gj — gj), where B,;y = dg.Pg.\g=g* , 
9i = {g^^}- The stability coefficients 9i, defined as minus the eigenvalues of By, provide an 
important characteristic of the fixed point. In particular, eigcndircctions with a positive (negative) 
real part 9 arc UV-attractive (UV-repulsive) for trajectories close to the fixed point. For the 
remainder of this subsection, we will set d = 4. 

Inspection of (2.37) immediately gives the Gaussian fixed point (GFP) 

g* = 0, X*=0, rjl, = rj* = 0. (3.1) 

This fixed point corresponds to the free theory, and constitutes a saddle point in the 5-A-plane. 
It has one attractive and one repulsive eigendirection with stability coefficients given by the 
canonical mass dimensions of G and A, respectively. 

The numerical analysis of the ghost-improved /3-functions also reveals a unique NGFP. Its 
position and properties are shown in the first two lines of Table 1. It is situated at ^r* > 0, A* > 
and UV-attractivc in both g, A. Substituting its position into r]c determines the ghost anomalous 
dimension rj* = —1.8.^ For comparison, the third line of Table 1 displays the properties of the 
NGFP obtained within the standard Einstein-Hilbert truncation without ghost-improvement. We 
observe that the actual numerical values of the physical product g*X* and the stability coefficients 

''When including tiie marginal in the set of coupling constants, rj* is also the stability coefficient associated 
with the new (UV-irrelevant) eigendirection. Here we refrain from adopting this viewpoint, however, since the 
running of tJc is completely determined by gk, Afc so that is, most likely, not an essential coupling. 
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Truncation 


A* 


9* 


g*\* 




Re(^) 


Im(6l) 


cutoff 


EH + ghost 
EH + ghost 


0.135 
0.260 


0.859 
0.355 


0.116 
0.092 


-1.774 
-1.846 


1.935 
2.070 


2.012 
2.439 


opt 
exp (s = 1) 


EH 


0.193 


0.707 


0.136 




1.475 


3.043 


opt 



Table 1: Properties of the NGFP arising from the ghost-improved /3-functions (2.37). The first 
two hnes show the position, the universal product g*X*, the ghost anomalous dimension r/*, 
and the stability coefficients of the fixed point obtained with the optimized cutoff (A. 17) and 
exponential cutoff (A. 19) with s = 1, respectively. For comparison, the third line displays the 
characteristics of the NGFP found in the standard Einstein-Hilbert trunaction [27] . 



are shifted by approximately 30%. This is in the typical range for the cutoff-scheme dependence 
observed in [27]. Most remarkably, both the standard and the ghost-improved Einstein-Hilbert 
truncation give rise to the same fixed point structure. This is highly non-trivial, as the new 
contributions to the gravitational /3-functions are of the same order of magnitude as the other, 
already known, terms. We interpret this result as a striking confirmation of the gravitational 
fixed point structure disclosed by the standard Einstein-Hilbert truncation [17, 18, 15, 16, 27]. 

The main virtue of the ghost-improvement becomes apparent, when investigating the stabil- 
ity of the physical quantities 5* A*, ^, rj* with respect to the variation of the IR cutoff R^- To 
illustrate this point, we resort to the exponential cutoff (A. 19), and determine the properties 
of the NGFP for varying shape parameter s. Fig. 1 shows the resulting cutoff-scheme depen- 
dence of the physical quantities for the standard (dashed red line) and ghost-improved (solid 
blue line) computation, respectively. Remarkably, the ghost-improvement reduces the unphysical 
cutoff-scheme dependence by factors 1.4, 2.0 and 3.2 for the product g*X* ■, Re^ and Im^, respec- 
tively. The scheme-dependence of rj* can only be determined in the ghost-improved computation. 
Here the variation is approximately 2%. Since the cutoff-scheme dependence provides the pro- 
totypical probe for judging the reliability of a truncation, these results clearly indicate that the 
ghost-improvement significantly improves the quality of the Einstein-Hilbert approximation. In 
particular the fixed point properties given in Table 1 should be more robust as the ones obtained 
without ghost-improvement. 

3.2 The ghost-improved phase portrait 

After analyzing the fixed point structure, we determine the phase portrait resulting from the 
ghost-improved /3- functions. We start by investigating the gravitational RG flow on the g- 
A— plane, before studying the behavior of the anomalous dimensions along some typical sample 

trajectories. 

The phase portrait resulting from the numerical integration of the /3-functions (2.37) is de- 
picted in Fig. 2. We first observe that = is a fixed line, which cannot be crossed by the 
flow. For 5fe > the flow is dominated by the interplay of the NGFP and the GFP. In this 
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Figure 1: Stability analysis for universal quantities in the Einstein-Hilbert truncation (dashed red 
line), previously obtained in [15], and upon including the ghost wave-function renormalization 
(solid blue line). The ghost-improvement significantly decreases the cutoff-scheme dependence. 



regime, the UV behavior of the RG trajectories is controlled by the NGFP, which acts as a UV 
attractor. Following the RG flow from this fixed point towards the IR, the RG trajectories un- 
dergo a crossover from the NGFP to the "classical regime" dominated by the GFP. Depending 
on whether the trajectory turns to the left (Type la), right (Type Ilia) or hits the GFP (Type 
Ila), the classical theory has a negative, positive or zero cosmological constant. The trajectories 
with positive cosmological constants, however, cannot be continued to k = 0, but terminate at 
a finite value of k when reaching the boundary of the phase space. The latter is indicated by 
the green line, which constitutes a singularity in the /3-functions at finite g, A. Observe that the 
ghost-improved phase portrait, Fig. 2, is in complete qualitative agreement with the standard 
Einstein-Hilbert truncation [15, 27, 8]. 

It is now illustrative to pick one sample trajectory for each of the classes discussed above and 
study the A;-dependence of the anomalous dimensions tjn and rjc along the flow. This is shown in 
Fig. 3. In the UV, for large values of t = ln(A;/A;o), the anomalous dimensions are determined by 
the NGFP, so that ry^ = — 2 and rj* = —1.77. Lowering t and approaching the IR, the anomalous 
dimensions undergo a crossover towards the classical theory with ?7jv ~ 0, ?7c ~ 0. The steep 
increase at the end of the Type Ilia trajectory is caused by the singularity of the /3-functions 
(green line in Figure 2), and heralds the termination of the trajectory at a finite value k. The UV- 
limit of the ratio rjc/rj^ is governed by the NGFP and takes the value ri*/rj'^ = 0.89. Following 
the flow towards the IR, the ratio undergoes a crossover and asymptotes to 0, the cutoff-scheme 
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-0, 4 



-0,2 



Figure 2: The RG flow of g^, \k obtained from the numerical integration of the ghost-improved 
/3-functions (2.37) using the optimized cutoff. The green line indicates a boundary of the coupling- 
constant space where the /3-functions diverge. The phase portrait is in complete agreement with 
the one obtained from the standard Einstein-Hilbert truncation [15, 27, 8]. 

dependent value 54$3'7(24$2 - 

or a finite value at the termination point, for trajectories 
of Type la, Type Ila, and Type Ilia, respectively. 

It is worthwhile to have a closer look at the singularity causing the termination of the Type Ilia 
trajectories. Fig. 3 shows, that at this point in coupling constant space, the anomalous dimensions 
r/jv and 77c diverge. This can be traced back to the vanishing of the denominators in (2.39). Here, 
the ghost-improvement has, however, a very non-trivial effect: while the denominator arising 
from the standard Einstein-Hilbert truncation has a term linear in g only, the ghost-improvement 
adds an additional piece quadratic in g. This may provide an elegant mechanism for lifting this 
singularity by shifting the zeros of the denominator to complex values g. However, for the Bi, Ci 
given in Appendix B, this mechanism is not realized, and may requires a further improvement of 
the truncation before becoming operational. 

3.3 The NGFP in space-times with extra-dimensions 

The /3-functions (2.38) and (2.39) are continuous in the space-time dimension. In the sequel, we 
will exploit this feature and analyze the resulting gravitational fixed point structure for general d. 
Based on the standard Einstein-Hilbert truncation, a similar analysis has already been carried out 
in [15, 19, 21]. Motivated by the recent interest in asymptotically safe TeV-scale gravity models 
[43, 44, 45, 46], which hinge on the existence of the NGFP in the presence of extra-dimensions, 
it is worthwhile to complement the previous results by including the ghost-improvements. 

Remarkably, the d-dimensional fixed point structure is strikingly similar to the one obtained 
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Figure 3: RG flow of tlie anomalous dimensions r^Tv (upper left), ijc (upper right), and their ratio 
Vc/vn (lower right) along the sample RG-trajectories of Type la (orange). Type Ila (black), and 
Type Ilia (blue) given in the lower left diagram. 

in four dimensions. The GFP (3.1) exists for all d. Furthermore, there is a unique generalization 
of the NGFP for all dimensions 3 < d < 25 considered here. Its properties are summarized in 
Fig. 4. The fixed point is situated at positive 5* > 0, A* > and UV attractive in both g and 
A. Below d < 24 its critical exponents are given by a complex pair Re0 it ilm^. For d > 24 the 
imaginary part of the critical exponents vanishes and we have two real stability coefficients, which 
are still UV attractive. All these results are in perfect agreement with earlier findings based on 
the standard Einstein-Hilbert truncation. An interesting feature arises in the ratio ^*/??7v, with 
r/^ = 2 — d, shown in the top right diagram of Fig. 4. This ratio is peaked at d = 4, where 
it reaches almost unity. For both d > 4 and d < 4 the ratio decreases. It is a rather curious 
observation that in d = 4 the graviton and ghost propagators have the same anomalous dimension 
in the UV, highlighting this particular value of the space-time dimension.^ While it is clearly 
desirable to get a better understanding of this result, the underlying analysis is beyond the scope 
of the present paper. 

4 Discussion and Conclusions 

In this paper we have analyzed the gravitational renormalization group flow arising from the 
ghost-improved Einstein-Hilbert truncation. Besides a scale-dependent Newton's constant and 

related argument, concluding that d = 4 is special, is based on the spectral dimension of space-time computed 
from the running of Newton's constant and the cosmological constant and has been put forward in [47] . 
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Figure 4: Properties of the NGFP for general space-time dimension d, obtained with the optimized 
cutoff (A. 17). There is a unique NGFP for all 3 < d < 25 considered. 

cosmological constant, this setup also takes into account the quantum effects captured by the 
non-trivial wave-function renormalization in the ghost sector. Our main result is the occurrence 
of a unique non-Gaussian fixed point of the gravitational /3-functions. It is located at positive 
Newton's constant and cosmological constant and UV attractive in both g and A. Its properties 
are strikingly similar to the non-Gaussian fixed point found in the standard Einstein-Hilbert 
truncation [15, 16, 17, 18, 19, 20, 21]. This finding further substantiates the asymptotic safety 
scenario of quantum gravity [1].^ The fixed point also persists for higher-dimensional space- 
times, solidifying the foundations for asymptotically safe TeV-scale gravity models [43, 44, 45, 46], 
possibly accessible at the LHC . 

The main virtue of the ghost-improvement is a significant decrease of the unphysical cutoff- 
scheme dependence in physical observables as, e.g., the critical exponents of the fixed point or the 
universal product g*X*. Since this scheme dependence constitutes a standard test for the quality 
of a truncation, we conclude that the ghost-improvement significantly increases the robustness 
of the emerging physical picture. In particular, the characteristics of the fixed point displayed in 
Table 1 should provide a better approximation of the stability coefficients emerging from the full 
theory, as compared to the standard Einstein-Hilbert result. 

The ghost-improvement hinges on the new formula (2.39) for the ghost anomalous dimension 

^Qualitatively, our picture is also confirmed by the very recent results [41], which study the ghost-improved 
Einstein-Hilbert truncation employing a spectrally adjusted cutoff. The numerical variations observed in the two 
computations are within the typical range expected from the different cutoff-schemes [27]. 
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rjc, which is analytic in Newton's constant and completely determined by g and A. At the 
non-Gaussian fixed point, we obtain rj* = —1.77 for the optimized and ry* = —1.85 for the 
exponential cutoff. Notably, these numbers are very close to the anomalous dimension of Newton's 
constant 77^ = —2. We believe that this is not accidental, but a consequence of space-time 
becoming effectively two-dimensional at short distances, where the physics is controlled by the 
non-Gaussian fixed point [47, 48]. The essence of our argument is the observation that the 
propagator of a field with anomalous dimension r] is proportional to p~'^~^^. Therefore, in the 
vicinity of the UV fixed point, both the graviton and ghost propagator behave approximately 
as p^^, which translates into a logarithmic correlator in position space. In this case the space- 
time seen by both fields is effectively two-dimensional, suggesting rj* = —2 in the full theory. 
A similar spontaneous dimensional reduction has also been observed within the framework of 
Causal Dynamical Triangulations [49], and, by now, also in a variety of other quantum gravity 
approaches [50] . It would be very interesting to explore, if this is a first hint pointing at a deeper 
connection between these seemingly unrelated fields. 
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A Heat-kernel techniques 

In this appendix we collect the essential formulas for evaluating the operator traces (2.28) and 
(2.29). We start with the traces containing only Laplacians before turning to operator insertions 
including the background ghost fields in Subsection A. 2. Their relation to the threshold-functions 
employed in the main text is established in Subsection A. 3. 

A.l Early time heat-kernel expansion 

The differential operators entering into the traces (2.28) appear in the form of covariant Laplacians 
A = —D^ only, and may thus be evaluated using standard early time heat-kernel methods 
for second order differential operators (see [27] for more details). In this case the heat-kernel 
expansion takes the form 

TVi [e-^^] = (47rs)-'^/2 j (fx^ [tr^ ao + s tr^ 02 + ■ • • ] , (A.l) 
with coefficients oq = 1? and a2 = \R'^i and trj denoting a trace over the internal indices 

tr2Tl2T = 5(d + 2)(o!- 1), trill =d, trolo = 1 ■ (A.2) 
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Here and in the following the dots indicate interaction terms which do not contribute to the 
truncation. 

By employing a Laplace-transform, eq. (A.l) can be generalized to suitable operator- valued 
functions W{z), 

Tvi[W{A)] = {A7r)-^/^ J d^x^{Qa/2[W]tnao + Qd/2-i[W]ivia2 + ---} , (A.3) 

where 

POO -I pOG 

Q„[W]:= dss-^W{s) = —- dzz^'-^Wiz). (A.4) 
Jo ^ I'^j Jo 

Here W{s) is the Laplace transform of W{z) and we used a Mellin-transform to re-express the 

integral over W{s) in terms of the original function. Employing (A. 15) to recast the functionals 

Qn arising from (2.28) in terms of the dimensionless threshold functions (A.13), a straightforward 

computation leads to the right-hand-side of eq. (2.32). 



A. 2 Applying the Anselmi-Benini perturbation method 

While the traces Si can be evaluated using the standard heat-kernel techniques for Laplace op- 
erators, dealing with the Qi is more involved. Here, the main technical difficulty results from the 

fact that their arguments contain covariant derivatives which do not organizing themselves into 
minimal second-order differential operators. Applying the perturbative hcat-kernel technique 
developed by Anselmi and Benini [40], this obstruction can be circumvented by the following 
two-step procedure: In the first step, wc use commutator relations to collect all non-minimal 
operators in a single operator insertion. The second step then uses the heat-kernel formula at 
non-coincident points to evaluate the operator traces including these "perturbations" . 

To implement the first step, we start from the arguments in (2.29) and express the propagator 
sandwiched between the A's (or Q's) by its Laplace transform. Then, all the Laplacians reside in 
an exponential, that can easily be commuted with the A's (or Q's) 

[O, e-'"^] = ~sc-'^[0, A] + 0(s2) . (A.5) 

Undoing the Laplace-transform leads to the desired form of the trace argument, including only 
one operator of non-Laplace-type. The prefactor s multiplying the commutator terms thereby 
translates into a derivative of the original propagator. At the level of the full trace argument, 
(A.5) implies 

W{A)Q'"'a = Q^'aWiA) - 2W'{A) (D^Q^^'a) + • • • , (A.6) 

and similarly for the other terms. Since our sole interest is in isolating the contribution propor- 
tional to the ghost kinetic term, all terms containing more than two covariant derivatives acting 
on and or involve curvature tensors do not contribute to the truncation and may therefore 
be dropped. 

After these manipulations, the traces Qi assume the generic form 

POO 

Tri[OW{A)]= dsW{s){x\Oe-''^\x) . (A.7) 
Jo 
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Here all Laplacians have been collected in PF(A) and O symbolizes a (matrix- valued) insertion 
build from two covariant derivatives acting on the background ghosts together with two (or four) 
covariant derivatives acting to the right. (See below for an explicit example.) The trace (A.7) is 
evaluated by inserting a complete set of states between the two operators 

{x\Oe-'^\x) = {x\0\x'){x'\e-'^\x) = J d'^x^iii [O H{s,x,x')\^,^^ . (A.8) 



, oo 

cr(x,x ) 



Here 

H{s, X, x') := {x'\e-'^\x) = (Atts)-'^/'^ e"^^ ^ s''A2n{x, x') , (A.9) 

n=0 

is the expansion of the heat-function at non-coincident points, with A2n{x,x') being the off- 
diagonal heat-kernel coefficients and a{x, x') denoting half the squared geodesic distance between 
x and x' , satisfying 

u{x,x) = Q, ^a'/^a-fj, = a , (T.^y{x,x) = g^^{x) , a.(„^...„^)(a;, x) = , n > 3 . (A.IO) 

Notably, A2n{x,x) = a2„, so that all but Aq{x,x) give rise to curvature terms. Furthermore, 
covariant derivatives acting on A2n{x,x') either vanish or produce higher derivative curvature 
terms once the coincidence limit is taken [51, 40]. Thus, when applying O on the function 
if(s, X, x'), the only terms contributing to the truncation are those where the covariant derivatives 
act pairwise on a{x^x'). 

To illustrate the procedure, let us consider the following explicit example from the A°^Aa 
product in Qq 

Oo = [\{D^e)D^][-\{D,c^)D^] = -j,{D^c"){D,Ca)D^D'' + • • • . (A.ll) 

Here we dropped all terms with more than two covariant derivatives acting on the background 
ghosts. Substituting Oq into (A.7) gives 

Tr[W{A)Oo] = —-^-^rdss-'''^-^W{s) f d^x^{D,c'^){D^^c^) 



where we used (A. 4) in the second step. Applying this algorithm systematically to all terms in 
(2.29) and subsequently using the relation (A.16) then leads to (2.33). 

A. 3 Threshold functions and cutoff scheme 

In order to capture the cutoff-scheme dependence and highlight the structure of the /3-functions 
it is convenient to express the functionals (5n[l^] in terms of dimensionless threshold functions. 
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The structures appearing in the Qi motivate the definitions 



1 r 
1 r 

W)Jo 

1 r 
1 



dz z 



dz z 



dz z 



dz z 



[z + i?(o) {z) + wY {z + R(o) (z)) ^ ' 

n-l 

{z + Ri'^){z) + wY {z + R(o){z)y ' 
{z + i?(o) {z) + wY {z + (z)) ^ ' 



(A.13) 



n-l 



r(n) 7o " {z + Rio) (z) + wY {z + (z)) « ' 

where the prime denotes the derivative of the shape function defined below (2.19), with 

respect to its argument. Observe that the threshold functions with p = do not depend on w. 
These definitions naturally generalize the ones used in [12], which are recovered as special cases 

$^'V) = ^^H, <'^H = $^(o), $^'°h = $^;h, 4>0'^h = |.^(o). (A.14) 

These threshold functions are closely related to the Q-functionals (A. 4). For the traces (2.28) 
we have 

dtiZlRk) 



Qn 



ZliPk + w)p\ 



2^2(n-p+l) $P,0^^/^2^ _ 1 ^ ^P^%w/k'^) 



(A.15) 



with the index I = N,c such that the relation holds for both the gravitational and ghost wave- 
function renormalization. This relation generalizes to the Q-functionals occurring in (2.29) 



Qn[m=2Zlk 



c ;„2ri-4 



ir/e^i'2(-2Afe: 



(A.16) 



Qn[f2 



■■2Zl.k 



I , 2n-6 



$2-2(_2A,) + $f (-2A,) - i,7,(*f (-2Afc) + $f (-2Afe)) 



In combination with the heat-kernel formulas (A. 3) and (A. 7), these equations allow us to find 
the /3-functions for Z^,Z^ and Afe, (2.32) and (2.35). 

We close this section with a remark on the shape functions employed in the numerical studies 
of Section 3. Unless stated otherwise, all results are obtained with the optimized cutoff [52] 



i?f°P*(z) = (l-z)G(l-z), 



(A.17) 



In this case the integrals appearing in (A.13) can be carried out analytically 

1 1 . 1 1 



r{n + I) {I + w)P ' ^ r(n + 2) (1 + w;)P ' 



(A.18) 



r(n + 1) (1 + u')P ' ' r(n + 2) (1 +1^)^ ■ 

Here, the threshold functions degenerate such that they become independent of the index q. When 
analyzing the cutoff-scheme dependence of our truncation, we also employ a one-parameter family 
of smooth exponential shape functions 



^^■^^■^(z;^) 



sz 



exp(sz) 



(A.19) 
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The continuous shape parameter s allows to smoothly vary the implementation of the IR cutoff. 
In contrast to the optimized cutoff, the integrals in the threshold functions cannot be carried 
out analytically for this class of cutoffs. Thus we have to resort to numerical integration when 
evaluating the threshold functions. 

B Functions determining the anomalous dimensions 

In this appendix, we give the definitions of the functions Bi{\) and Ci(A), completing the con- 
struction of the anomalous dimensions for the graviton and ghost fields (2.39). In the following, 
all threshold functions arc evaluated at the argument w = — 2A, which we then suppress for 
notational simplicity. The Bi are exactly the same as those obtained in [12] 

B,{X) = 1 (4vr)i-'^/2 [did + l)$;;,°_i - 6did - 1)$^;, - 4d^'^l,_, - 24$;;,) , 
B2{X) = - i(47r)i-<^/2 [did + + 6did - l)*^;,) • 

The quantum corrections from the wave-function renormalization of the ghosts are encoded in 
Ci(A) = (47r)^-'^/^ (2Cgr$^)Vi - 4rf(^S+2 + ^^72+2)) ' 

CiX) = I (4;r)^-'^/2 (2^4;^^, + 12l>2;,y , 
C4(A) = - (47r)i-<^/2 {4C^^l%^, + 4Cgh$if2+i) ' 
with the coefficients Cgr and Cgh defined in (2.36). 
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